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Abstract 

We present the g-Stirling's formula using the g-product determined by Tsalhs entropy as nonex- 
tensive generahzation of the usual Stirling's formula. The numerical computations and the proof 
are also shown. Moreover, we derive and prove some inequalities and equalities concerning the 
g-Stirling's formula. Applying the g-Stirling's formula, the mathematical foundations of Tsallis 
statistics can be obtained as similarly as those of Boltzmann-Gibbs statistics. 
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I. INTRODUCTION 



Tsallis statistics derived from Tsallis entropy [1] [2] provides many successful descriptions 
of the physical systems with power-law behavior, which cannot be systematically recovered 
in the usual Boltzmann-Gibbs statistical mechanics [3]. The appearance of such a new 
physics stimulates us to finding its corresponding new mathematical structure behind it [4] . 
On the way to the goal, we obtain the q- Stirling's formula as nonextensive generalization of 
the usual Stirling's formula. 

The g-Stirling's formula is derived from the q-product determined by Tsallis entropy. The 
g-product is first introduced by Borges in [5] for satisfying the following equations: 

lnq{x ^qU) ^lUqX + lnqy, (l) 

expg (x) ®g expq (y) = exp^ (x + y), (2) 

where In^ x and exp^ {x) are the q-logarithm function[3]: 

— 1 

luq X := — — — (x > 0, g e M+) (3) 
and its inverse function, the q- exponential function[3]: 

( [l + (l-q)x]^ iil + (l-q)x>0, 
I otherwise 

The requirements (1) and (2) lead us to the definition of ®q between two positive numbers 

[x^-'i + - 1]^ , if X > 0, 1/ > 0, + y^-'i - 1 > 

(5) 

otherwise 

which is called the q-product in [5]. The g-product recovers the usual product: 

lim {x ®q y) = xy. (6) 

The fundamental properties of the g-product ®q are almost the same as the usual product, 
but the distributive law does not hold in general. 

a{x®qy) ^ ax®qy (a e M+, x, e M+) (7) 

In our previous paper [6], a Tsallis distribution is derived from applying the g- product to 
the likelihood function in maximum likelihood principle (MLP for short) instead of the usual 




product. This first successful application of the g'-product to the MLP is briefly reviewed as 
follows. 

Consider the same setting as Gauss' law of error: we get n observed values: 

Xi,X2, • • • e M (8) 

as results of n measurements for some observation. Each observed value Xi {i — 1, ■ • • ,n) is 
each value of identically distributed random variable Xi {i — 1, ■ ■ ■ ,n), respectively. There 
exist a true value x satisfying the additive relation: 

= x + Cj (i = 1, • • • ,n) , (9) 

where each of Cj is an error in each observation of a true value x. Thus, for each Xi, 
there exists a random variable Ei such that Xi — x + Ei {i — 1, ■ • • ,n). Every Ei has 
the same probability density function / which is differentiable, because ^^i, • • • ,X„ are 
identically distributed random variables, (i.e., Ei,--- ,En are also identically distributed 
random variables.) Let Lg (6) be a function of a variable 9, defined by 

Lq (6) = Lq (Xi, Xs, • • • , X„; 6*) := / (Xi - 9) ®q f {X2 -9)®q---®qf{Xn-9). (10) 

If the function Lq{xi,X2, - ■ ■ ,Xn;9) oi 9 for any fixed Xi,X2, - ■ ■ ,Xn takes the maximum at 

9^9* + + + 
n 

then the probability density function / must be a Tsallis distribution: 

expg {—Pqe^) 

^ ^ Je^Pqi-Pqe^)de ^^^^ 

where /3q is a ^'-dependent positive constant. This result recovers Gauss' law of error when 
0' — > 1. See [6] for the proof. 

In addition to the above successful application of the gr-product, we show one more im- 
portant property of the g-product in section II. This property convinces us of the undoubted 
validity of the g-product in TsaUis statistics. Along the lines of the above successful applica- 
tion and the important properties of the g-product, we derive the g-Stirling's formula with 
the proof and numerical computations in section III. 
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II. REPRESENTATION OF THE g-EXPONENTIAL FUNCTION BY MEANS OF 
THE g-PRODUCT 



The base of the natural logarithm e is well-known to be defined as 

exp (1) := lim ( 1 + . (13) 

n-+oo y nj 

Using this definition, the exponential function exp (x) is expressed by 

exp (x) = lim f 1 + -) . (14) 

n^oo \ n/ 

Replacing the power on the right side of (14) by the n times of the g-product (g)^ : 

X®« ■.= X(^q---(»qX, (15) 
n times 

exp (x) is obtained. In other words, lim (1 + -) '' coincides with exp (x) . 

exp (x) — hm 1 H — (16) 

expg ix) in (4) was originally introduced in [7] from Tsallis entropy and its maximization. 

1 + - ) ''is computed using the definition of the g-product only. 
Therefore, this coincidence is an evidence revealing the conclusive validity of the g-product 
in TsaUis statistics as similarly as the law of error in Tsallis statistics. Therefore, the q- 
Stirling's formula using the g-product is naturally defined and computed in the following 
section. The proof of (16) is given in the appendix. 

III. g-STIRLING'S FORMULA 
A. gf-factorial nlq 

The g-factorial n\q for n e N and g > is defined by 

n\q := 1 (8)g • • • (8)g n (17) 
where ®q is the g-product (5). Using the definition (5), n\q is explicitly expressed by 



1-9 

n\q= ^ fc^-"^ - (n - 1) . (18) 
.fe=i 



This imphes 



n 



1-q 



(19) 



In the present paper, we derive the approximation of In^ {n\q) instead of that of n\g. In fact, 
when q > 1, there exists infinite numbers of n e N such that ^^^=1 k^~'^ — (n — 1) < 0. For 
example, 

n 

k^-" - (n - 1) < when q = 1.2 and n > 6. (20) 



fe=i 



When X^^^i A;^~* — {n — 1) < 0, n\q takes complex numbers in general except q — 1 + 
2^ (m e N). On the other hand, when J2k=i - (n - 1) > for a given n e N and 
g > 0, we easily obtain n\q from In^ {n\g). Therefore, the approximation of In^ {n\g) is called 
"g-Stirling's formula" in the present paper. If an approximation of lug {n\q) is not needed, 
the above explicit form (19) should be directly used for its computation. However, in order 
to clarify the correspondence between the studies q = 1 and q ^ 1^ the approximation of 
In^ {n\g) is useful. In fact, using the present g-Stirling's formula, we obtain the surprising 
mathematical structure in TsaUis statistics [8]. 
We obtain the ^'-Stirling's formula as follows: 



In^ {n\g) 



where 



(n + I) In n + (-n) + 6'„,i + (1 - ^i) 
^ -Inn- I + ^,,,2-^2 

_1_ ( ^ g 



n 



2-q 



n I 1 \ 1 — 1 



l-q 



2-q 



^ ) + On,q + ( fZ^ - 



lim en,q = 0, < On,! < 



a-<5i 



12n' 



if g = l, 
if g = 2, 

if g>0 and q 1,2, 

(21) 



(22) 



It is easily verified that the g-Stirling's formula coincides with the usual Stirling's formula 
when g — > 1. 

linilng (n!g) = Inn! (23) 

Other properties of 9n,q and Sq such as explicit forms and approximations are still unclear. 
But numerical results of (21) exhibit good approximations of the true value (19), which is 
shown in the next subsection. 
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B. Numerical computation of the g-Stirling's formula 



By the property lim 6'„^q = in (22), 6n,q is set to be for the computation. The 
expUcit form of 6q is not clear yet, but 6q does not depend on n. Thus, dq is set to be 
Sq — 5i — 1 — In V27r ~ 0.081 which is already known in the usual Stirhng's formula. Then 
the g-Stirling's formula for the numerical computations in this subsection is given by 



hiq (nlq) 



(n + I) In n+(-n) + (1-0.081) if q ^ 1, 

-- Inn- i- 0.081 if q ^ 2, 

+ i) ^ + (-^) + - 0.081) if q>0 and g^l,2, 



n 



(24) 

According to the following numerical results, these assumptions on 9n,q — and Sq — 0.081 
are found to be justified when n ^ oo. We compute the true value (19) and its approximate 
value (24) for In^ (nlq). Moreover, the relative error and the absolute error are also computed. 

I true value — approximate value] 



relative error := , 

true value 

absolute error := jtrue value — approximate value | . 



(25) 
(26) 



g = 0.1 



n 


In, {nlq) (18) 


approximate value (24) 


relative error 


absolute error 


2 


9.6230 X 10-^ 


8.8677 X 10-1 


7.8480 X 10" 


-2 


7.5521 X 10" 


-2 


5 


9.2143 


9.1455 


7.4604 X 10- 


-3 


6.8742 X 10" 


-2 


10 


3.9708 X 10^ 


3.9644 X 10^ 


1.6117 X 10" 


-3 


6.3995 X 10- 


-2 


15 


9.0005 X 10^ 


8.9943 X 10^ 


6.8181 X 10" 


-4 


6.1366 X 10" 


-2 


20 


1.5933 X 10^ 


1.5927 X 10^ 


3.7384 X 10" 


-4 


5.9563 X 10" 


-2 


i 






i 




i 




oo 

















(27) 
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71 


niq \n.q) yi-O) 


"nriT'n'viTTij^'l'p vj^lnp i i 

ClL/L/l W-AJ.±±lCl/uC veil Lie- 


relative error 


absolute error 


9 




7 71 1 9 V 1 n~l 


u.yiou A iu 


-2 


c; 7Q1 1 V 1 ri" 

0. ( O i i A iU 


-2 


o 


7fi47 




5 2937 X 10" 


-3 


3 5810 X 10" 


-2 


10 


2.4937 X 10^ 


2.4912 X 10^ 


9 9893 X 1 0' 


-4 


9 AQin X 10" 


-2 




5 0938 X 10^ 


5 0918 X 10^ 


3.9413 X 10- 


-4 


2.0076 X 10- 


-2 


20 


8.3332 X 10^ 


8.3315 X 10^ 


2.0633 X 10- 


-4 


1.7194 X 10- 


-2 


i 






i 




i 




oo 

















7i 


In, (,;!,) (18) 


a])])r()xiiiia.te value (21) 


rela.ti\'e error 


al)S()lute error 


2 


6.9339 X 10-1 


6.5208 X 10-1 


5.9570 X 10- 


-2 


4.1305 X 10- 


-2 


5 


4.7906 


4.7740 


3.4666 X 10- 


-3 


1.6607 X 10" 


-2 


10 


1.5118 X 10^ 


1.5110 X IQi 


5.4802 X 10- 


-4 


8.2851 X 10- 


-3 


15 


2.7930 X 10^ 


2.7924 X 10^ 


1.9711 X 10- 


-4 


5.5052 X 10- 


-3 


20 


4.2387 X 10^ 


4.2383 X 10^ 


9.7063 X 10- 


-5 


4.1142 X 10" 


-3 


i 






i 




i 




oo 

















n 


In (n!) 


approximate value (24) 


relative error 


absolute error 


2 


6.9315 X 10-1 


6.5187 X 10-1 


5.9553 X 10- 


-2 


4.1279 X 10- 


-2 


5 


4.7875 


4.7709 


3.4639 X 10" 


-3 


1.6583 X 10- 


-2 


10 


1.5104 X lOi 


1.5096 X lOi 


5.4746 X 10- 


-4 


8.2691 X 10- 


-3 


15 


2.7899 X lOi 


2.7894 X lOi 


1.9690 X 10" 


-4 


5.4933 X 10- 


-3 


20 


4.2336 X lOi 


4.2332 X lOi 


9.6960 X 10- 


-5 


4.1049 X 10- 


-3 


i 






i 




i 




oo 

















q = 1.5 



q = 2.0 
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relative error 


absolute error 




At 


5 S57Q X 1 0~^ 


5 5504 X 1 0^^ 


c; 94SQ V 1 n~ 


-2 


Q 0747 X 1 n~ 


-2 


O 




S 5975 


2 5856 X 10" 


-3 


9 1442 X 10- 


-3 


10 


9.9580 


9.9537 


4 3609 X 1 0" 


-4 


4 3496 X 1 0" 


-3 




1 7172 X 10^ 


1 7169 X 10^ 

X.I X \J i_/ /\ X w 


1.8303 X 10- 


-4 


3.1431 X 10- 


~3 


20 


2.4809 X 10^ 


2.4807 X 10^ 


1.0643 X 10- 


-4 


2.6406 X 10- 


-3 


i 






i 




i 




CX) 

















n 


In, (n!,) (18) 


approximate value (24) 


relative error 


absolute error 


2 


5.0000 X 10-1 


4.7585 X 10-1 


4.8294 X 10- 


-2 


2.4147 X 10- 


-2 


5 


2.7167 


2.7096 


2.6152 X 10- 


-3 


7.1046 X 10- 


~3 


10 


7.0710 


7.0664 


6.5292 X 10- 


-4 


4.6168 X 10- 


-3 


15 


1.1682 X 10^ 


1.1678 X 10^ 


3.5564 X 10- 


-4 


4.1545 X 10- 


-3 


20 


1.6402 X 10^ 


1.6398 X 10^ 


2.4342 X 10- 


-4 


3.9926 X 10- 


~3 


i 






i 




i 




oo 

















(31) 



(32) 



The numerical results (27)~(32) and Fig.l show that the relative error of the g-Stirling's 



0.08 
0.07 
0.06 



■3 0.05 

J 0.04 

2 0.03 

^ 0.02 

IS 

U 0.01 



10 



15 



20 



FIG. 1: relative error of the g— StirUng formula with respect to n when q = 0.1, 0.2, • • • , 2.0 
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formula (24) in any q > goes to zero when n — > oo. Therefore, the g-Stirhng's formula 
(21) can provide a good approximate value for each In^ (n!g) in any g > 0. 

C. Derivation of the g-Stirling's formula 

This subsection proves the g-Stirhng's formula (21) when g > (g 7^ 1). The case g = 1 
of the g-Stirling's formula (21) is the famous usual Stirling's formula whose proof is found in 

many references such as [9], so that wc prove the g-Stirling's formula in case g > (g 7^ 1). 
For simplicity, we assume g 7^ 1 in the subsections C and D. 

Note that, when g = 0, n!^ is directly and explicitly derived from (18) as follows: 

n\g=o = Y.k-{n-l)^-n''--n+l. (33) 

k=i 

Thus, the g-Stirling's formula when g = is not needed. When g < 0, In^ a; is a convex 
function with respect to x, so that Sn,q defined by (37) below does not converge when n — > 00 
in general. Moreover, Tsallis entropy when g < lacks the property of concavity which 
plays essential roles in the natural generalization of Boltzmann-Gibbs statistical mechanics. 
Therefore, we prove the g-Stirling's formula (21) when g > 0. 

Consider the definite integral J" In^ xdx when g > (g 7^ 1) . The curve in Fig. 2 describes 
a concave function with respect to x: 

y = lnqX= ^ _ ^ ■ (34) 

For each natural number A; e {2, 3, • • • } on the x-axis, there exist a rectangle with width 
1 and height In^ A; (See Fig. 2). Here, in order to avoid confusions, we call a rectangle at 
/c G {2, 3, ■ ■ ■ } on the x-axis by "/c-rectangle" . The center of the bottom of each /c-rectangle 
is at e {2, 3, • • • } on the x-axis. For each /c-rectangle, we define ak,q and Pk,q for e N by 

ak,q '■= the area of the portion surrounded by the curve, the top side of the /c-rectangle 
(x-axis when k — 1) and the left side of the {k + 1) -rectangle, (35) 

Pk,q '■— the area of the portion surrounded by the curve, the left side 

of the {k + 1) -rectangle and the top side of the {k + 1) -rectangle. (36) 

Then, we define the sum of areas of all gaps between the curve and the rectangles by 

9 



y = ln„x 




FIG. 2: approximation of In^ xdx 



n— 1 n— 1 71—1 

^Jijij • ^ ^ (^k,q 

k=l k=l k=l 



Pk,g = {ak,q - Pk,q) ■ 



(37) 



Using 6n,q, the area of the portion surrounded by the x-axis, the straight hne x = n and the 
curve y = lugX is computed as 



luq 2 + In, 3 H h luq (n - 1) + ^ In, n + 6n,q 



luq xdx 



" /xi-<? - 1 

1-g 



Inn + n — 1 if <7 = 2 
^■r.-ln,n-^ if 



2-9 



Therefore, we obtain 

lug = In, (1 (g). 



n) = lUg 2 + In, 3 + ■ • ■ + In, n 
1 



In, 2 + In, 3 H h In, (n - 1) + - In, n + 6n,q 



lnn + n-l) + i(l-i)-V2 if q = 2 



- In, n - dn,q 



^ n ■ ln„ n - ^ ] + i In, n - if q^2 

if g = 2 



2n 



^ + i)ln,n+ + if q^2 



(38) 
(39) 

(40) 

(41) 

(42) 
(43) 
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liig X is a concave function of x, so that the magnitude relation between ak,q and /3k,q is 
shown in Fig. 3. From Fig. 3, we find 

tangent line 



k,q at X = k + 1 





k+l,q 



y = ln„x 



FIG. 3: magnitude relation between a^^q and Pk,q 



> f3l,q > Ol2,q > ■■ ■ > > (^n,q > 



(44) 



This means that /3i,g, a2,g, ■ ■ ■ , an,q, f^n,q, ■ ■ ■ " is a monotone decreasing sequence with 

lower bound 0, so that we construct a sequence of terms of an alternating series such that 

~Pl,qi Ct2,q, " /^2,q, " " " , Otn,q, —l3n,q, ' ' ' ■ (45) 

Thus, 6n^q defined by (37) is found to be a finite partial sum of terms of the above sequence 
(45). Therefore, 

oo 

lim Sn,q = {ak,q - (3k,g) (46) 

n— »oo ' 

fe=l 

is an alternating series. Note that in general an alternating series is of the form 
YliLi i~^y~^^ ('^i > 0)- It is known that for an alternating series with property (44) Sn,q 
converges as n — > cxo for any fixed q > 0. Thus, we define 6q and 6'„ g as 



6q := lim 6„^g, 

ra— »oo 



^n,q • ^q ^n.q- 



(47) 
(48) 



Substituting 5„ g = Sq — 6n,q into (43) yields 



lUg {n\q) 



n - ^ - \nn - \ - {82 - Onfl) 

n I 1 
2-q 



if q = 2 



(49) 
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where 5q — lim Sn,g implies that 

n— >oo ' 



hm 9n,q = hm {6q - Sn,q) = 0. (50) 
n— >oo n—*oo 

Therefore, when n is sufficiently large, In^ (n!q) is approximated as 



■ n— J — Inn — ^ — S2 ifg — ^ 

The remained problem is to determine dq which does not depend on n. This is still an open 
problem, but the exact value of Sq for actual applications is not needed. The assumption 
Sq — 61 — 1 — In V2Tr gives sufficiently useful approximations of the g-Stirling's formula (21) 
as shown in numerical computations in the previous subsection. 



D. Inequalities and equalities concerning the g-Stirling's formula 

This subsection presents the following inequalities (52)~(53) and equalities (54)~(55) 
concerning the g-Stirling's formula. In some applications, these formulas are more useful 
than the g-Stirling's formula (21). 

n . n 1 , 1 1 , / , \ 
lUg n - h - lUg n + <hiq [n\q) 



2-q " 2-q 2 " 2-q 



Tl 72 1 1 

<2^^"^^"2^ + 2^^^"+2^ 



(52) 

15 11 
^-Inn- — - - < lnq(n!<,) < n-lnn- — - - (if ? = 2) , 

(53) 

where n e N and g > 0. These inequalities or the g-Stirling's formula (21) is reduced to the 
equalities: 

In, (77!,) = ^ In, 77- 23^ + (In, 77) (if q ^ 2) , (54) 

In, (77!,) = 77 - In 77 + O (1) (if q = 2). (55) 

These formulas are proved in the following way. 
In Fig. 4, the sum of gray triangular areas is 

^ In, 2 + i (In, 3 - In, 2) + ^ (In, 4 - In, 3) + • • • + ^ (In, 77 - In, (77 - 1)) = ^ In, 77. (56) 

12 



y 



In (n - 1) 



y = ln,j X 




FIG. 4: J" liiq xdx and the sum of the triangular areas 



On the other hand, 



Ing {n\q) = Ing 2 + In,, 3 H h In^ tt, 



(57) 



imphes that In^ {n\q) represents the sum of the rectangular areas in Fig.4. Thus, from Fig. 4 
we can read that the difference between lug {n\q) and /" In^ xdx is less than the sum of gray 
triangular areas. 

(58) 



In, (ri!g) 



luq xdx < - luq n 



Therefore, we obtain an upper bound of In^ {n\g) as follows 

r 1 

lUg (nig) < / luq xdx + -InqU 
Ji 2 

(-Inn + n- 1) + i (1 - i) if q = 2 
^-^-n-ln^n-ffi) + iln,n if q^2 



ri — Inn 



1 1 



if q = 2 



^ln,n-^ + iln,n + ^ if q^2 



(59) 
(60) 

(61) 



2-q 



where we used (39). 



We derive an lower bound of \nq{n\q). 
Ik~^ -'^^^ G {2, ■ ■ ■ ,n — 1}) is written by means of /9fc-i,g and a^,, 



In Fig. 5, the integral 



luq xdx = \ngk- f3k-i,q + otk,q (A; G {2, ■ ■ ■ , n - 1}) 



(62) 
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y 




FIG. 5: upper bound of /" liig xdx 



Using the relation shown in Fig. 3: 



ak,q is replaced by (3k-i,q in (62), so that we get 

In, xdx <\nqk- (3k-i,q + [3k-i,q = liig k. 

Summing up the both sides of (64) for k G {2, ■ ■ ■ , n — 1}, 

_ 1 

2 

luq xdx < lug 2 + ■ ■ ■ + luq (n — 1) . 



a'l q in Fig.5 is a triangle surrounded by the tangent line of y = In^x at x 



and a: = |. Thus, 



On the other hand, 



In, xdx < a[ , 



ln„ xdx < - ln„ n. 

-I 2 



The sum of each side of (65), (66) and (67) yields 

r 1 1 

/ lug xdx < - + In, 2 + ■ ■ ■ + lUg (n — 1) + - lUg n 
Ji 8 2 



In, {n\g) + i - ^ In, n. 



(63) 
(64) 

(65) 

1, the X-axis 
(66) 
(67) 

(68) 
(69) 
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Thus, 



luq xdx + Q In, ~ < i'^U) ■ i^^) 
Substituting (39) into J^lngxdx on the left side of (70) yields 

ln„ (nlo) > < / \ (71) 

(^•n-Hn-f5j) + (|Hn-|) if q^2 



n — Inn— ^ — I ifq' = 2 



_^ln^^-^ + iln,n + ^-i if 
Prom (61) and (72), we obtain 

n , n 1 1 1 , / , N 

H^- 7^ — :+ 0^*^+ 7^ — - - o < K (n!g) 



(72) 



2-g " 2-q 2 ^ 2-q 

n , n 1 , 1 /.r / r.\ 

< hian 1 — lUonH (if o#2), 

2-q * 2-q 2 " 2-q ^ 



(73) 
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n — In n < Ino in\a) < n — In n (if o = 2) 

2n 8 ^' 2n 2 ^ ^ ' 



(74) 



This result also shows 



H(ng--— ln,n--— + 0(ln,n) (if g ^ 2) , (75) 
lug {n\g) = n - In n + O (1) (if g - 2) . (76) 



IV. CONCLUSION 

We present the g-Stirling's formula in Tsallis statistics with the proof and some numerical 
results. Moreover, some inequalities and equalities concerning the g-Stirling's formula are 
shown. The present g-Stirling's formula can be applied to many scientific fields in Tsallis 
statistics as similarly as the usual Stirling's formula. In particular, the g-Stirling's formula 
plays important roles in Tsallis statistics when we want to clarify the difference between 
the studies q — 1 and g 7^ 1. Some typical and important applications of the g-Stirling's 
formula in Tsallis statistics will be presented in our forthcoming paper [8] , which reveals the 
surprising mathematical structure behind Tsallis statistics. 
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APPENDIX A: THE PROOF OF (16) 



In this appendix, we prove that exp (x) is derived from hm (1 + -) ^ using the q- 
product only. 



Using the g- product (5), (l + ^) 



1 + - 
n. 



1-9 



is exphcitly expanded as 



n I 1 + - 
n 



(n-1). 



The power series expansion for the right side of (Al) about the point ^ = is 



n 



X 1 



X\' 
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n 



Taking n ^ oo on the both sides yields 



lim 

n— ►oo 



n I 1 + - 
n 



(n-1) 



l + {l-q)x. 



Therefore, if 1 + (1 - g) x > 0, then 



lim 

n— >oo 



lim 

n— *oo 



n(l + -j -(n-1) 
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[1 + (1 - q) x] 1-" 
exp^ (x) . 
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